Structure and Melting of Two-Species Charged Clusters in a Parabolic Trap 
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We consider a system of charged particles interacting with an unscreened Coulomb repulsion in a 
two-dimensional parabolic confining trap. The static charge on a portion of the particles is twice 
as large as the charge on the remaining particles. The particles separate into a shell structure with 
those of greater charge situated farther from the center of the trap. As we vary the ratio of the 
number of particles of the two species, we find that for certain configurations, the symmetry of 
the arrangement of the inner cluster of singly-charged particles matches the symmetry of the outer 
ring of doubly-charged particles. These matching configurations have a higher melting temperature 
and a higher thermal threshold for intershell rotation between the species than the nonmatching 
configurations. 

PACS numbers: 82.70.Dd, 64.60.Cn, 83.10.Mj 
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Clusters of repulsive particles in confined traps have 
attracted considerable recent attention due to their ap- 
plicability to a wide variety of systems. For example, 
two-dimensional (2D) clusters can represent electrons in 
quantum dots [1] or on the surface of liquid helium [2] , 
vortices in superfluids [4], colloidal particles in circu- 
lar traps [5], confined ferromagnetic particles [6], and 
charged dust particles in plasma traps [7]. The 2D 
charged clusters also resemble the problem of charge dis- 
tribution studied by Thomson in his "plum-pudding" 
model of the atom [8] . 

When confined to a parabolic trap, charged particles 
form a structure of concentric rings, with the inner par- 
ticles forming a distorted triangular lattice resembling a 
defected Wigner crystal, and the outer rings taking on a 
more circular shape that conforms to the radial symme- 
try of the trap [9,10]. Among the possible charge config- 
urations are several "magic" arrangements in which the 
number of particles is such that the shells form with only 
a few symmetrically distributed dislocations, and so have 
a reduced total energy compared to what is predicted 
based on a semi-empirical approximation [11]. Bubeck 
et al. [5] have observed that certain colloidal clusters 
confined by a circular hard-wall trap exhibit two-stage 
melting, in which intershell rotation between the outer 
two shells occurs at temperatures below the temperature 
at which particles are able to jump between the shells. 

Several explanations of this two-stage melting phe- 
nomenon have been proposed [9,12-14], all of which focus 
on the intershell rotation which occurs prior to the ex- 
change of particles between shells. Most plausible among 
these is the theory that the rotation is due to an incom- 
mensuration between the shapes of the potentials created 
by the adjacent shells. For this intershell rotation to 
occur, the inner shell configuration must be sufficiently 
stable to have a melting temperature higher than the 
threshold for intershell rotation. 



In our simulation, we extend the confined charge sys- 
tem to include particles with two distinct values of 
charge. We find that the two species separate into shells, 
with those of greater charge located farther from the cen- 
ter of the trap. This occurs due to the fact that the 
pinning force couples only to the position of the parti- 
cle, while the interparticle repulsion is a function of both 
position and charge. Thus particles of greater charge 
are pushed farther up the walls of the parabolic trap. If 
two-stage melting and intershell rotation is caused by an 
incommensuration between the potentials formed by the 
particles, then we should be able to predict the occur- 
rence of intershell rotation based on the ratio of particles 
in the outer shell of singly-charged particles to the num- 
ber of doubly-charged particles. 

We consider a system of N s + Nd charged particles 
interacting via an unscreened l/r Coulomb repulsion, 
where N s is the number of single-charge particles with 
q s = 1 and is the number of doubly-charged par- 
ticles with qd — 2. The particles are free to move in 
two-dimensions but are confined by a parabolic trap cen- 
tered at the origin and increasing radially as r 2 . The 
dimensionlcss Hamiltonian [10] for this system is: 
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where qi {fi) is the charge (position) of particle i and 
we fix the strength of the parabolic trap to A = 10. 
Charged colloidal particles in a confinement potential, 
such as those employed by Wei et al. [3], bear the clos- 
est resemblance to the Coulomb interaction used in this 
study. 

Using a molecular dynamics (MD) simulation method, 
we initialize the system at high temperature, simulated 
by random Langevin kicks, and then slowly anneal it to 
a T = ground state configuration. We checked the 
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FIG. 1. Ground-state colloid configurations for: (a-d) 
N d = 5, (e-h) N d = 6, and (i-1) N d = 7, with N s = (a) 
6, (b) 8, (c) 11, (d) 31, (e) 7, (f) 12, (g) 19, (h) 30, (i) 8, (j) 
19, (k) 22, (1) 24. Small dots are the singly charged particles 
and large dots are the doubly charged particles. 

accuracy of this method by reproducing the ground state 
configurations for single-species particle clusters found 
by Kong et al. [15], albeit with different confinement 
strength. Once we have obtained the ground states of 
the two-species clusters, we slowly increase the tempera- 
ture and observe the melting of the system. 

Fig. 1 shows several examples of the ground-state 
configurations we obtained for = 5, 6, and 7, respec- 
tively. For Nd=5, configurations with N s = 6, 11, and 31 
have five-fold symmetry matching the number of doubly- 
charged particles in the outer shell. N s = 8 approxi- 
mates a six-fold symmetric arrangement by substituting 
a singly-charged particle to fill a potential well in the 
outermost shell. For Nd = 6, configurations with N s = 7 
and N s = 19 have six-fold symmetry. Configurations 
with N s = 12 and N s = 27 have three-fold symmetry 
and configuration N s = 30 has two-fold symmetry, which 
also are commensurate with the arrangement of particles 
in the outer shell. For iV^ = 7, N s = 8 and 22 are the 
only seven-fold symmetric inner particle configurations. 
N s = 23 and 24 are nearly symmetric, however, with 
seven-fold symmetry in the outer two shells of the inner 
particle configuration. In general we find highly ordered 
structures when the smaller particles form a commen- 
surate structure with the outer particles. All the con- 
figurations are shown at [16]. The general expressions 
N s = kNd and N s — kNd + 1, with k a small integer, 
predict some configurations that have rotational symme- 
try through an angle 2n/N d , such as N d — 5, N s = 6, 
but do not work for all values of k due to the fact that 
the outer particles may distort from a uniformly spaced 
arrangement in order to better accommodate the inner 
particles. Although in this work we only present results 




T 

FIG. 2. Plot of <j s vs. temperature for N d — 6, N a — 7, 
averaged over five realizations. 



for qd = 2q s , we have considered other charge ratios, 
and find that the more highly charged particles always 
move to the outside of the trap, and the same general 
commensurate-incommensurate phases occur. We have 
also performed simulations with larger Nd > 7; how- 
ever, we do not treat this case here since the outer par- 
ticles begin to form multiple rings which have their own 
commensurate-incommensurate transitions. In addition, 
we have run a set of simulations for N s < 15 with the 
same trap but with a 1/r 3 interparticle interaction poten- 
tial, which corresponds to the interaction between mag- 
netic colloids. We found that, in general, the ground 
states and also the dynamical properties were qualita- 
tively identical to those in the case of Coulomb repulsion 
for small N s (N s < 15). 

In order to confirm whether the highly ordered com- 
mensurate phases are more stable we consider the melting 
of the two-species system. We determine the tempera- 
ture T e of the first exchange of particles between shells 
by measuring intershell exchange of the singly charged 
particles, 

1 Ns 

^ = ^E^(*)-^(°)i ( 2 ) 

i—l 

This gives the mean radial distance of the inner singly- 
charged particles from their initial T = positions. a s 
shows a marked increase at T e when the inner particles 
begin to jump between shells, as illustrated in Fig. 2 for 
a system with N s — 7 and Nd = 6. This configuration 
is highly stable, as shown in Fig. 1(e), with one central 
particle surrounded by an inner hexagonal shell of six 
singly-charged particles and an outer hexagonal shell of 
six doubly-charged particles. At T = 3.55, a a jumps 
when one of the first-shell particles exchanges with the 
central particle, as shown in the trajectories of Fig. 3. 

A similar measure, ad, tracks the doubly-charged par- 
ticles; however, as they never formed more than one shell 
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(a) (b) 
FIG. 3. Particle trajectories (lines) for Nd = 6, N s = 7 at 
temperatures (a) T — 3.5 with no intershell exchange and (b) 
T — 3.6 with intershell exchange. 



for the parameters considered here, there were no in- 
tershell exchanges and ad only increased linearly with 
temperature. We observe exchanges of a singly-charged 
particle with a doubly-charged particle only at very high 
temperatures r > T e , and these exchanges occur only 
for highly asymmetric and disordered configurations. 

To track the onset temperature T r of intershell rotation 
between the two species when it occurs, we use a second 
measure, Ag: 

1 N a 1 N d 

A * = \w - ^°)) - w 5> j(<) _ ^ (0))l ' 

s i=l j=l 

(3) 

which gives the difference between the mean angular dis- 
placements of the two species from their initial configura- 
tions. Ag increases when the shells slip past each other, 
but is insensitive to coherent rotation of the two species. 
Ag becomes meaningless if the particles do not maintain 
the same orientation with respect to the other particles of 
the same species, so it can detect intershell rotation only 
when this occurs before the onset of intershell exchange. 
Ag also detects relative slip between shells of the same 
particle species, which does not interest us here. Such 
same-species shell slips were generally limited to the er- 
ratic rotation of two particles at the center of the config- 
uration, and produced a sufficiently small change in Ag 
to be distinguished easily from a genuine rotation rela- 
tive to the other species. For configurations with a single 
central particle located roughly at the origin, the angu- 
lar displacement of that particle was excluded from Ag 
in order to reduce noise. 

An example of intershell rotation measured by Ag is 
shown in Fig. 4, where we plot Ag vs temperature for 
a configuration of N s — 7 and Nd — 7. In this case, the 
configurations of each of the two species are highly stable 
independently, with the singly-charged particles forming 
a hexagonal ring around a central particle. However, 
there is an incoherence between the hexagonal ring of 
singly-charged particles and the seven-particle ring of 




FIG. 4. Plot of A g vs. temperature for N d = 7, N s = 7, 
averaged over five realizations. 




(a) (b) 

FIG. 5. Particle trajectories over a short time interval for 
Nd = 7, N a — 7 at (a) T — 0.6 with no intershell rotation and 
(b) T = 0.7 with intershell rotation. 



doubly charged particles which surrounds it. A substan- 
tial increase in Ag occurs at T rs 0.65, corresponding to 
a slipping between the outer shell and the inner cluster 
of particles. The trajectory plot, Fig. 5, also indicates 
this slipping, as the trails of the inner shell particles be- 
come larger than those of the outer shell particles at this 
temperature. Since the onset of intershell rotation at T r 
is more gradual than the intershell exchange transition, 
we define T r to occur when Ag exceeds a threshold value 
of Ag = jj-, corresponding to an angular displacement 
of half of the angular distance between neighboring par- 
ticles in the outer shell of singly-charged particles. 

The matching of the inner particle symmetry with 
the outer shell symmetry produces an elevated melting 
threshold for commensurate configurations. The melting 
temperature T m is taken to be the lower of T e or T r . In 
Fig. 6(b) we plot T m for a range of N s at fixed Nd = 6. 
Figs. 6(a) and (c) show T m for Nd — 5 and Nd = 7, re- 
spectively. For each configuration, we averaged o~ s and 
Ag over five realizations to reduce error. In Fig. 6(b), 
T m for Nd — 6 has a peak at N s — 19, another perfect 
triangular arrangement. We also find an 
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FIG. 6. Melting temperature T m vs N s for (a) N d 
N d = 6, and (c) Nd = 7. 
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elevated melting temperature for Nd — 12, which has 
threefold symmetry. As seen in Fig. 6(a), configuration 
Nd = 5 and N s = 11 also exhibits a higher than average 
melting temperature. This is expected since the config- 
uration N s — 11 has five- fold symmetry, commensurate 
with the number of outer particles. Fig. 6(c) shows high 
T m values for N d = 7 at N s = 23 and N s = 24, which 
is expected due to the matching seven-fold symmetry of 
the outer three rings. The inner part of the configuration 
Nd = 7, N s = 19 forms the same highly stable perfectly 
hexagonal arrangement of Nd — 6, N s = 19, and thus 
has a high threshold for intcrshcll exchange between the 
singly-charged inner particles. However, the incommen- 
suration with the 7 particles on the outside yields a low 
threshold for intershell rotation. 

There is an obvious shift in the peaks of Fig. 6(a- 
c) as the number of outer-shell doubly-charged particles 
changes. Most notably, the highest melting temperature 
for each Nd occurs for N s = Nd + 1, when the singly- 
charged particles form a central particle surrounded by 
a shell that is commensurate with the outer particles. 
As expected, the highest melting temperature occurs for 
Nd = 6 and N s = 7, as this arrangement matches the 
triangular lattice of the Wigner crystal. However, melt- 
ing temperatures nearly as high occur for Nd — 5 and 
N s — 5 and 6. We believe that the smaller size of the 
crystal in this case exaggerates the effect of the commen- 
suration between rings, despite the dissimilarity with the 
bulk lattice configuration. 

In conclusion we have investigated the structure and 
melting of two species of charged particles in a parabolic 
trap. The more highly charged particles form an outer 
ring. Highly ordered clusters occur when the structure of 
the central cluster of singly charged particles is commen- 
surate with the outer ring of doubly charged particles. 
We observe variations in the melting temperatures of the 
two species clusters, with elevated melting thresholds for 



intershell rotation and intershell particle exchange when 
a commensuration occurs between the symmetry of the 
inner cluster and the outer ring. 

This work was supported by the U.S. Department of 
Energy under Contract No. W-7405-ENG-36 and by the 
U.S. DoE, Office of Science, under Contract No. W-31- 
109-ENG-38. BJ and JD were supported by NSF-NIRT 
award DMR02-10519 and the Alfred P. Sloan Founda- 
tion. 



[1] M.A. Reed and W.P. Kirk, Nanostructure Physics and 
Fabrication (Academic Press: Boston, 1989). 

[2] P. Leiderer, W. Ebner, and V.B. Shilkin, Surf. Sci. 113, 
105 (1987). 

[3] Q.-H. Wei, C. Bechinger, D. Rudhardt, and P. Leiderer, 

Phys. Rev. Lett. 81, 2606 (1998). 
[4] Y. Kondo, J.S. Korhonen, M. Krusius, V.V. Dmitriev, 

E.V. Thuneberg, and G.E. Volovik, Phys. Rev. Lett. 68, 

3331 (1992). 

[5] R. Bubeck, C. Bechinger, S. Neser, and P. Leiderer, Phys. 

Rev. Lett. 82, 3364 (1999). 
[6] M. Golosovsky, Y. Saado, and D. Davidov, Phys. Rev. E 

65, 061405 (2002). 
[7] W-T. Juan, Z-H. Huang, J-W. Hsu, Y-J. Lai, and Lin I, 

Phys. Rev. E 58, R6947 (1998). 
[8] J.J. Thomson, Philos. Mag. 7, 237 (1904); B. Partoens 

and F.M. Peeters, J. Phys.: Condens. Matter 9, 5383 

(1997). 

[9] V.M. Bedanov and F.M. Peeters, Phys. Rev. B 49, 2667 

(1994) . 

[10] V.A. Schweigert and F.M. Peeters, Phys. Rev. B 51, 7700 

(1995) . 

[11] A. A. Koulakov and B.I. Shklovskii, Phys. Rev. B 57, 
2352 (1998). 

[12] V.A. Schweigert and F.M. Peeters, Phys. Rev. B 51, 7700 
(1995). 

[13] I.V. Schweigert, V.A. Schweigert, and F.M. Peeters, 

Phys. Rev. Lett. 84, 4381 (2000). 
[14] Y-J. Lai and Lin I, Phys. Rev. E 64, 015601 (2001). 
[15] M. Kong, B. Partoens, and F.M. Peeters, Phys. Rev. E 

65, 046602 (2002). 
[16] All configurations can be seen at 

http: / /www.t 12.lanl.gov /home/olson/Colloid2S .html. 



4 



